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Magnetization Process of the S = 1 and 1/2 Uniform and Distorted 
Kagome Heisenberg Antiferromagnets 
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The magnetization process of the 5 = 1 and 1 /2 kagome Heisenberg antiferromagnet is studied 
by means of the numerical exact diagonahzation method. It is found that the magnetization 
curve at zero temperature has a plateau at 1/3 of the full magnetization. In the presence of 
\/3 X \/3 lattice distortion, this plateau is enhanced and eventually the ferrimagnetic state is 
realized. There also appear the minor plateaux above the main plateau. The physical origin of 
these phenomena is discussed. 
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§1. Introduction 

The kagome Heisenberg antiferromagnet (KHAF) has 
been extensively studied theoretically and experimen- 
tally because of the interest in the interplay of the strong 
quantum fluctuation and the highly frustrated nature of 
the lattice structure. So far, most of the attempts have 
been focused on the ground state and low Mna^excita- 
tions of the uniform KHAF. In both S = l/'Mm^ and 
S* = 1E|) cases, it is known that the ground state is a spin 
singlet state and the magnetic excitation has a finite en- 
ergy gap. In the 5 = 1/2 case, there are a number of 
singlet excitations below the first triplet excitation pos- 
siblly down to zero energy in the thermodynamic limit D-' 
On the other hand, the singlet excitations also have fi- 
nite energy gaps in the 5* = 1 case. The present author 
proposed the hexagonal singlet salid (HSS) picture for 
the ground state S ^ I KHAF.B 

In the present work, we first investigate the magneti- 
zation process of this model. As examples of the frus- 
trated quantum magnets, the magnetization process 
the triangular Heisenberg magnet is widely studied.BB 
The magnetization plateau at one-third of the full mag- 
netization is well established in the Ising-like classical 
Heisenberg model as atypical example of the 'classical' 
magnetization plateau.EP Although this plateau can be 
well understood in terms of the classical spin configura- 
tion, it has been observed even in the isotropic Heisen- 
berg modeLwith S = 1/2 which has strong quantum 
fluctuation.B' On the other hand, the honeycomb lattice 
S* — 1/2 Heisenberg antiferromagnet has no plateau-El* 
In this context, it is interesting to investigate the case of 
kagome lattice which has the triangles and hexagons at 
the same time. __ _ 

As a real material, Wada and coworkersEil'lld'Ej' have 



This material, however, undergoes a structual transfor- 
mation at-i28.7K to the distorted phase with -^/S x -s/S 
structure.tj^ Motivated by this observation, we further 
investigate the effect of such lattice distortion to the 
ground state and magnetization process in both S* = 1/2 
and S* = 1 cases. We find that the plateau is enhanced 
by the lattice distortion and eventually the ferrimagnetic 
ground state is realized. 

This paper is organized as follows: In §2, the model 
Hamiltonian is presented. The numerical results for the 
magnetization curves are presented in §3. The physical 
pictures of the plateau state and ferrimagnetic state are 
also explained. The last section is devoted to a summary 
and discussion. 

§2. Model Hamiltonian 

The Hamiltonian of our model is given by, 



N 



(2.1) 



where H is the magnetic field, g the electronic g-factor, 
fj,B the Bohr magneton and Si is the spin operator. In 
the following, we take the unit g/ie = 1. The summation 
S<i j> is taken over all the nearest neighbour pairs of 
sites of the kagome lattice. 

We also consider the V3 x V3 distorted version of this 
model given by. 



7i = Ti-A + "Hb + "He 



<i.,j>Ea 



N 
i=l 



(2.2) 



(2.3) 



' i- u 1 • c n^^nAAATATm^ 1- 1 where V^, represent the summation over the bonds 
mvestigated the magnetic behavior of m-MPYNN-BF4 which i u '-^ ^ a td u u- i. 



can be regarded as the 5 = 1 kagome antiferromagnet. 



around the type-a (a 
depicted in Fig. 1. 



A, B or C) hexagons, which are 
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Fig. 1. \/3 X V3 distorted kagome lattice. 
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§3. Magnetization Curves 

The Hamiltonians of the finite size clusters shown in 
Fig. 2 are diagonahzed numericaUy to obtain the zero 
temperature magnetization curves. The periodic bound- 
ary condition is assumed. 

3.1 Uniform case 

For S — 1/2, we use the clusters with N ~ 12, 15, 
18, 21, 24 and 27 to draw the full magentization curve. 
The clusters with iV = 30 and 33 are used only to calcu- 
late the high field region of the full magnetization. For 
5 = 1, we use the clusters with N — 12, 15 and 18 
to draw the full magentization curve. The cluster with 
TV = 21 is used at higher field region. The magnetization 
curves of the S" = 1/2 and S = 1 cases are shown in Fig. 
3. Apparently, there exists a magnetization plateau with 
magnetization AI = Ms/3 where Ms(= NS) is the satu- 
ration magnetization. The end points of the plateaux are 
plotted against the system size in Fig. 4. It is evident 
that the size dependence of the width of the plateau is 
weak in the S = 1 case while it is rather strong in the 
5=1/2 case. This would be due to the localized nature 
of the HSS state reahzed in 5 = 1 KHAF.B However, 
even for the 5=1/2 case, the width does not strongly 
depend on the system size for N > 21. Therefore we 
may expect that the 1/3-plateau survives in the thermo- 
dynamic limit both for^ =1/2 and 5=1. 
Recently, ZhitomirskyH has pointed out that the mag- 




Fig. 3. The magnetization curves of (a) 3 = 1/2 {N = 
12,15,18,21,24,27,30 and 33) and (b) S = 1 KHAF (A^ = 
12, 15, 18 and 21). Here and in the following figures the mag- 
netic field is normalized by J/gfiB- 



H (b) 



Fig. 2. Clusters used for numerical diagonalization. 
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Fig. 4. The system size dependence of the position of the plateau 
of (a) S = 1/2 and (b) 5 = 1 KHAF. The lower (upper) end of 
the plateau is denoted by Hci{Hc2)- 



netization plateau at Ms/3 in the classical KHAF is sta- 
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bilized by the order-from-disorder mechanism due to ther- 
mal fluctuation. Similar phenomenon is also pointed out 
in the classical triangulaiU'' andpfrustrated square lat- 
tice Heisenberg antiferromagnet.ta-' In these cases, the 
plateau state consists of the collinear "ft J, structure on 
the triangular plaquettes, which is selected by the large 
entropy gain due to the soft modes associated with this 
structure. It should be noted that the classical tti states 
themselves are still infinitely degenerate due to infinite 
number of tiling patterns of ffj structure on the trian- 
gular plaquettes. It is known that the quantum fludjua,- 
tion also plays a similar role in the triangular lattice.El'E2P 
Therefore this mechanism should also work in quantum 
KHAF. In the quantum case, the ground state is stabi- 
lized because of the low zero point energy of the above 
mentioned soft spin wave modes. In analogy with the 
classical case, the plateau state in the quantum KHAF 
can be regarded as the quantum superposition of the "f "f | 
tiling structures. 

From this point of view, the absence of the plateau for 
the honeycomb lattice is easily understood. Although 
the honeycomb lattice, as well as the kagome lattice, can 
be regarded as a depleted triangular lattice, all triangles 
are depleted in the honeycomb lattice so that no plateau 
remains. 

3.2 Distorted kagome lattice 

There are various ways to introduce distortion in the 
present model. To be specific, we concentrate on the 
following three cases: (A) Ja = J{1 + S), Jb = J,Jc = 
J(l - 6) , (B) Ja = J, Jb = J(l - 6), Jc = J(l - 5) 
and (C) Ja = J, Jb = J(l — <^), Jc = 0- The numerical 
calculation is performed with = 18 and 27 for 5=1/2 
case and with A^ = 18 for 5 = 1 case. 
3.2.1 Case (A); Ja = J(l + 6), Jb = J, Jc = ^(1 - S) 

The magnetization curves are shown in Fig. 5 for S — 
1/2 and Fig. 6 for S — 1, respectively. It is evident 
that the plateau region is enhanced as S increases and 
eventually the ferrimagnetic state is realized for large 
enough S. 

The enhancement of the plateau and the appearance 
of the ferrimagnetic state for large S is understood as 
follows. For large 5, the 6 spins around each A-hexagon 
approximately form a 6-spin singlet state. The number of 
remaining spins is A^/3. The effective coupling between 
these spins are weak so that the magnetization of these 
spins easily saturate. To magnetize further, large energy 
to break up the singlet on A-hexagon is required. Thus, 
the plateau is enhanced by distortion. Especially, the ef- 
fective interaction between the 'alive' spins mediated by 
the second strongest bond B is ferromagnetic. Therefore, 
if Ja ^ Jb ^ Jc, the 'alive' spins form ferrimagnetic 
state with magnetization M — Mg/S. 

It should be remarked that this structure of the plateau 
state is consistent with the picture of the plateau state 
of the undistorted KHAF which is the quantum super- 
position of the tti tiling structures. In the plateau state 
of the strongly distorted KHAF, one spin remains 'alive' 
for each triangle and all of them are up-spins. The re- 
maining two spins on the same tringle are coupled via 



JA-bond and antifcrromagnctically correlated. This is 
a kind of superposition state of tti states although the 
position of one of the up-spin is strongly limited. There- 
fore, the plateau state and the ferrimagnetic state in the 
distorted KHAF are continuously connected with that in 
the undistorted KHAF. 
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Fig. 5. The magnetization curves of 5 = 1/2 distorted KHAF 
with (a) S = 0.2 and (b) 0.5 for Af = 18 and 27 for case (A). 




Fig. 6. The magnetization curves of 5 = 1 distorted KHAF with 
(a) S = 0.2 and (b) 0.5 for Af = 18 for case (A). 



3.2.2 Case (B).- Ja = J, Jb = J(l - S), Jc = J(l - S) 

In this case, the ground state remains singlet as shown 
in Fig. 7 and Fig. 8. The physical picture of the ground 
state in this case is discussed in detail in ref. ^ which 
is summarized as follows. In the strong distortion limit 
Ja ^ Jb — Jc, the strongest effective interaction among 
'alive' spins is the next nearest neighbour antiferromag- 
netic interaction according to the perturbational calcula- 
tion in rcf. 
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If other weaker effective interactions are 
neglected, the whole lattice of 'alive' spins are decom- 
posed into three equivalent sublattices of kagome type. 
These 'alive' spins again form singlet states on large 
hexagons which are three times larger than the original 
kagome lattice aad the ground state of the whole system 
is again singlet .llZP 

The slight backbending of the magnetization curve is 
observed above the plateau for the N — 27 cluster with 
S = 1/2,(5 = 0.2 (Fig. 7(a)). This implies the meta- 
magnetic behavior, if it remains in the thermodynamic 
limit. This phenomenon, however, takes place only in a 
narrow paramater range around 5 = 0.2 of the N ~ 27 
cluster. Within the present calculation, therefore, it is 
not clear if this is the artifact of the finite size effect or 
truely thermodynamic effect. This problem is left for 
future studies. 
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Fig. 7. The magnetization curves of 5 = 1/2 distorted KHAF 
with (a) <5 = 0.2 and (b) 0.8 for iV = 18 and 27 for case (B). 
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Fig. 8. The magnetization curves of 5 = 1 distorted KHAF with 
(a) 5 = 0.2 and (b) 0.8 for = 18 for case (B). 



3.2.3 Case (C);Ja = J,Jb = J(1 - S), Jc = 
In contrast to case (B), the ground state is always ferro- 
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Fig. 9. The magnetization curves of S = 1/2 distorted KHAF 
with (a) S = 0.2 and (b) 0.8 for A^ = 18 and 27 for case (C). 




Fig. 10. The magnetization curves of S = 1 distorted KHAF with 
(a) S = 0.2 and (b) 0.8 for Af = 18 for case (C). 



magnetic as shown in Fig. 9 and Fig. IQ—This is easily 
understood by the Lieb-Mattis theorem.EB' For Jc = 0, 



the whole lattice can be decomposed into two sublattices. 
One of the sublattice has 2N/3 sites and the other sub- 
lattice, which is coupled antiferromagnetically with the 
former, has sites. As a result, the ground state has 
total magnetization NS/3. This is also consistent with 
the argument in the strong distortion limit at the end 
of the last paragraph of case (A), because the effective 
interaction between the 'alive' spins are all fcrrimagnetic 
for Jc = 0. 

For distorted cases, there are minor plateaux at — 
(1/3 + 2m/9)A/s,m = 1,2 for S* 1/2 and at M, = 
(1/3 + m/9)Ms,rn ~ 1,..,5 for S* = 1 above the main 
plateau at AIz = Ms/3, although it is not clear if they 
remain in the thermodynamic limit. Of course, it is ob- 
vious that these plateaux appear for Jb, Jc ^ J A, be- 
cause these plateaux corresponds to the magnetization 
steps of isolated hexagons. Actually, in the undistorted 
case, such plateaux do not exist. Surprisingly, however, 
they also appear even if Jb and Jc are close to Ja as 
far as Jb = Jc- This is understood in the following 
way. If Jb — Jc, the effective field due to the polariza- 
tion of 'alive' spins are uniform on the six sites around 
the A-hexagon, so that it only cause the shift of effec- 
tive magnetic field and the positions of the plateaux are 
essentially those for the isolated hexagons. If the dif- 
ference between Jb and Jc is large, the effective field 
is nonuniform and the 6-spin eigenstates with different 
magnetization around the A-hexagon are strongly mixed 
up and the plateaux are smeared. In case (A), it should 
be noted that such plateaux are smeared for S = 1 but 
not for 5 = 1/2 for (5 = 0.2. Presumably, this is due to 
the stronger quantum effect in the S* = 1/2 case. In the 
5* = 1/2 case, the magnetization steps for the isolated 
hexagon is strongly quantized compared to the 5 = 1 
case so that the steps for the S* = 1/2 case are not easily 
smeared out. If the distortion is weaker, however, the 
minor plateaux are smeared out even in the S* = 1/2 
case. Therefore this difference is not an essential differ- 
ence between 5* = 1/2 and S = 1 cases but is rather a 
quantitative one. 

§4. Summary and Discussion 

The magnetization process of the spin- 1/2 and 1 KHAF 
is studied by means of the exact diagonalization method. 
The magnetization plateau at 1/3 of the full magnetiza- 
tion is found. The physical picture of the plateau state 
is explained as a quantum superposition of the ||| tiling 
states. It is shown that this plateau is enhanced by the 
X V3 lattice distortion and the fcrrimagnetic ground 
state is realized for large enough distortion in cases (A) 
and (C). The physical origin of the enhancement of the 
plateau and appearance of the fcrrimagnetic state is ex- 
plained from the argument in the strong coupling regime. 
It is also argued that the plateau state in the strong dis- 
tortion limit is continuously connected with that of the 
uniform case. The interpretation of the bahavior of the 
minor plateaux in the high field regime is also given. 

It should be interesting to observe the magnetization 
plateau experimentally. For to-MPYNN-BF4, the effec- 
tive exchange coupling between the S = l_£pins is given 
by J ~ 0.78K so that J/g^B 0.58r.BBB) Thus the 



Uniform and Distorted Kagome Heisenberg Antifcrromagncts 



5 



magnetic field at the plateau is easily feasible. The width 
of the plateau is around J. Therefore, the plateau would 
be clearly observed if the magnetization measurement at 
temperatures much below IK is carried out. 

The author thanks M. E. Zhitomirsky for attracting 
the author's attention to ref. 15, 16| and for useful com- 
ments to the earlier version of this work. The numerical 
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Supercomputer Center, Institute for Solid State Physics, 
the University of Tokyo and the HITAC SR8000 at the 
Information Processing Center, Saitama University. The 
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coded by H. Nishimori and KOBEPACK/1 coded by T. 
Tonegawa, M. Kaburagi and T. Nishino. This work is 
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from the Ministry of Education, Science, Sports and Cul- 
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